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ABSTRACT 

We present new, high resolution 1420 and 408 MHz continuum images and H I and ^^CO (J=l— 0) 
spectral line maps of the diffuse supernova remnant CTB104A (G93.7— 0.3). Analysis of the complex con- 
tinuum emission reveals no significant spectral index variations across the remnant. Three prominences 
around CTB104A are found to be related to the SNR, while one extension to the east is identified as an 
H II region associated with a background molecular shell. Small scale polarization and rotation measure 
(RM) structures are turbulent in nature, but we find a well-ordered RM gradient across the remnant, 
extending from southeast to northwest. This gradient does not agree with the direction of the global 
Galactic magnetic field, but does agree with a large-scale RM anomaly inferred from rotation measure 
data by Clegg et al. (1992). We show that the observed morphology of CTB104A is consistent with 
expansion in a uniform magnetic field, and this is supported by the observed RM distribution. By mod- 
eling the RM gradient with a simple compression model we have determined the magnetic field strength 
within the remnant as Bo w 2.3 /iC. We have identified signatures of the interaction of CTB104A with 
the surrounding neutral material, and determined its distance, from the kinematics of the H I structure 
encompassing the radio emission, as 1.5 kpc. We also observed clear breaks in the H I shell that corre- 
spond well to the positions of two of the prominences, indicating regions where hot gas is escaping from 
the interior of the SNR. 

Subject headings: HII regions — ISM: bubbles — ISM: individual (CTB104A) — magnetic fields — 
polarization — supernova remnants 
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1. INTRODUCTION 

CTB104A (G93.7-0.3) is a diffuse, thick-shelled super- 
nova remnant (SNR) with four "prominences" , extend- 
ing beyond the main body of the shell. Three of the 
prominences are polarized, as reported by Mantovani et al. 
(1991). The radio spectrum of CTB104A (S^ cx is 
relatively flat, typical of a mature shell- type SNR in the 
radiative expansion phase. 

In broad terms, CTB104A appears as a double limb- 
brightened shell. The lack of symmetry between the limbs 
may be indicative of different shock conditions across the 
remnant and hence a range of effective evolutionary ages. 
Kesteven & Caswell (1987) interpret CTB104A as a possi- 
ble barrel-shaped SNR with the shell of the remnant being 
physically distorted by the interstellar medium and elon- 
gated in the direction of the Galactic magnetic field. This 
is consistent with the observed morphology of CTB104A, 
but only for a particular choice of magnetic axis. Con- 
sequently, good quality polarization measurements can be 
used to examine this possibility. 

Apart from identifying the origin of the shell morphol- 
ogy, a complete understanding of CTB104A depends upon 
gaining information on the nature of the prominences ex- 
tending from the shell. With the exception of the promi- 
nence to the east (higher longitude) the prominences do 
not terminate at a sharp boundary, but rather appear 
to fade smoothly into the Galactic background emission. 
These prominences may be "hot blasts" escaping into re- 



gions where the density of the environment is presumably 
relatively low, as suggested by Landecker et al. (1985). A 
close examination of the surrounding neutral medium may 
reveal signatures relating to the origin of the prominences. 

To diagnose the dominant processes that shape 
CTB104A and to gain insight on the prominences extend- 
ing from the remnant, we have compiled high resolution 
radio continuum data at 1420 and 408 MHz, polarization 
measurements at 1420 MHz, H I spectral line data and 
HIRES processed IRAS images (Cao et al. 1997) from the 
Canadian Galactic Plane Survey (CGPS), as well as ^^CO 
(J=l— 0) emission line data at comparable resolution from 
the Five College Radio Astronomy Observatory (FCRAO). 
We examine the spectral index distribution with TT-plot 
analysis between 408 and 1420 MHz for signs of differ- 
ing shock conditions over the remnant's shell and also to- 
wards the prominences. Linear polarization measurements 
are used to study the overall role of the magnetic field 
on the expansion of the SNR and to search for spatially 
variable polarization that is indicative of regions of turbu- 
lence and material heated by localized shock compression. 
We identify signatures of the interaction of the SNR with 
its surroundings from the CGPS H I line data and use 
this additionally to constrain the distance to CTB104A 
as ^ 1.5 kpc. Finally, we use the FCRAO CO data and 
HIRES IRAS data to identify one of the prominences as a 
background H II region, located in the Perseus Arm, and 
thus not associated with the SNR. 
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2. OBSERVATIONS 

2.1. CGPS data 

CTB104A was observed at 1420 and 408 MHz in con- 
tinuum mode, including the Stokes Q, U and V parame- 
ters at 1420 MHz. Radio continuum and H I line obser- 
vations were carried out simultaneously with the Domin- 
ion Radio Astrophysical Observatory (DRAO) synthesis 
telescope (Landecker et al. 2000) as part of the CGPS. 
HIRES-processcd IRAS data (Cao et al. 1997) were also 
obtained from the CGPS data archive. The CGPS is de- 
scribed by Taylor et al. (2001) and a more detailed de- 
scription of the data processing routines can be found 
in Willis (1999). The angular resolution of the DRAO 
synthesis data varies slightly across the final maps as 
1' X 1' cosec((5)(1420/i^(GHz)). Single antenna data are in- 
corporated into the synthesis maps to assure accurate rep- 
resentation of all structures up to the largest scales. The 
low spatial frequency continuum data are obtained from 
the Effelsberg Survey at 1420 MHz (Reich et al. 1997) and 
at 408 MHz from Haslam et al. (1982). The low spatial 
frequency H I data are from the Low Resolution DRAO 
Survey of the CGPS region observed with the DRAO 26-m 
Telescope (Higgs & Tapping 2000). 

2.2. FOR AO CO data 

We have obtained spectral hne images of ^^CO (J=l— 0) 
emission in the area around CTB104A from the FCRAO 
14m telescope. Our mapped area covers Galactic longi- 
tudes ^ ~ 93° to ^ ~ 94?9 and Galactic latitudes 6 ~ -1° 
to & ~ -|-0?33, with a somewhat irregular bomidary. The 
CO data for £ < 94? 5 were obtained by Routledge and 
Moriarty-Schieven (unpublished) in November-December 
1994 using the QUARRY 15 pixel focal plane array (Er- 
ickson et al. 1992); to this we have added more data at 
e > 94?5, taken during March/April 2001, using the SE- 
QUOIA 16 pixel focal plane array (Erickson et al. 1999). 
Pointing measurements were carried out using Venus (1994 
data) and the SiO maser x Cygni (2001 data). Both sets 
of CO data were obtained by position-switching. 

The 1994 data have 80 MHz (~190 km s"^) total band- 
width, with the spectrometer centered on ulsr = —45 
km s"^ and the 2001 data have 40 MHz (~95 km s"^) to- 
tal bandwidth centered on vlsr = ^30 km s~^. The spa- 
tial sampling and velocity resolution (1.21 times the chan- 
nel spacing) for the 1994 and 2001 data sets respectively 
were 25", 0.98 km s'^ and 22", 0.25 km s'^. The spatial 
resolution (beam FWHM) of the data is 45". Both data 
sets were initially calibrated by the chopper wheel method 
(Kutner & Ulich 1981) and subsequently converted to the 
radiation temperature (T^) scale by correcting for forward 
scattering and spillover losses {rjfss = 0.7). Each data set 
was convolved to 2' spatial resolution and 0.98 km s~^ 
velocity resolution prior to combining them on the same 
^j^, ^^LSR grid as the H I data. The resulting sensitivity 
was 0.29 K and 0.08 K {la) for the 1994 and 2001 data 
sets respectively. 

2.3. Morphology of the SNR 

The radio continuum total intensity maps around 
CTB104A are shown in Figures 1 and 2, at 1420 and 

^ 20 - 40 degrees counter clock wise with respect to the GalcLctic plane. 



408 MHz respectively. The remnant has a roughly circu- 
lar shape consisting of three distinct fragments. A circle 
of 1?2 X 1?2 centered at ^ = 93?7 and b = -0?3, fitted 
by eye and shown in Figure 3, delineates this roughly cir- 
cular shape. The prominences, identified in Figure 3 by 
arrows, form an irregular boundary to the otherwise circu- 
lar remnant. The angular size of the SNR, taking all the 
prominences into account, is about 1?3 x 1?5. The north- 
ern part of the shell has the highest intensity at both 408 
and 1420 MHz maps. The prominence extending toward 
the south of the SNR, where the shell is incomplete, is a 
good candidate for an example of the "break out" phe- 
nomenon observed in mature SNRs (e.g. Cygnus Loop, 
VRO 42.05.01 and CTA 1). A similar claim could be made 
for the other prominences if projection effects are consid- 
ered. 

Disregarding the prominences, we examine the idea that 
the SNR's morphology could be interpreted as due to the 
expansion of the remnant in a regular magnetic field. In 
such a model (e.g. van der Laan 1962) no emission in 
total-intensity along the field direction is expected if the 
surrounding magnetic field is uniform on a scale larger 
than the SNR. The positions of the three shell fragments 
can support this model only for a particular choice of mag- 
netic field axis. This preferred axis, lying ~ 20 - 40 degrees 
north of west ^, partitions the remnant into two halves, 
one consisting of the north and east fragments and the 
other consisting of the southwest fragment alone. No other 
choice of magnetic field axis can validate this model. In 
Section 2.5 we show that the actual magnetic field axis 
lies at ^ 30° north of west. Thus, expansion in a uniform 
magnetic field is consistent with the observed morphology 
of the SNR. 

The diffuse filamentary structure with a bright core seen 
to the south of the SNR is the H II region S124. The dis- 
tance of SI 24 based on the distance of the exciting star 
from spectrometric data is d = 2.6±0.6 kpc (Felli & Harten 
1981, Brand & Blitz 1993), which excludes a physical re- 
lationship between SI 24 and the remnant since CTB104A 
is at a distance of ~1.5 kpc (see Section 3). 

2.4. Spectral analysis 

Spectral indices were computed using the conventional 
TT-plot method by plotting the brightness temperature 
values of the 408 MHz maps against those at 1420 MHz. 
The 1420 MHz data are convolved to the resolution of the 
408 MHz map. To avoid oversampling in the TT analy- 
sis the maps are regridded to 1'. In order to reveal any 
variation in the spectral index across the whole SNR, we 
have selected regions towards distinctly separable parts of 
the remnant and performed TT-plot analysis on the data 
extracted from those individual regions. Regions corre- 
sponding to the core of SI 24 and the extragalactic source 
4C50.55 were also included in this analysis for a reliability 
check on our measured indices. Figure 3 shows the selected 
regions marked on the 1420 MHz image. The regions were 
selected in a way to avoid unresolved sources, but this was 
not always possible. Since these sources are most likely 
extragalactic, we have checked the images for unresolved 
sources in the highest resolution image (1420 MHz) and 
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wc excluded a circular region with a diameter of 4' around 
these sources. Then we plotted brightness temperature 
values of both maps with respect to each other and fit- 
ted the distribution with a straight line. The slope of 
this line gives the temperature spectral index /3, where 
Tb ~ i^^ = We repeated this analysis by changing 

the independent variable and obtain another value for (3. 
We then adopted the difference of these two values of /3 
as the error in the spectral index. Spectral indices derived 
this way are not strongly affected by the presence of un- 
known backgrounds contributing to the emission bright- 
ness at each frequency. The presence of uniform back- 
grounds will modify only the unimportant constant term 
in the fit, leaving the derived slope unaffected, while vary- 
ing backgrounds will tend to increase the scatter. 

The results of the TT analysis are plotted in Figure 4. 
The analysis gives a = —0.1 for S124, which is typical for 
an optically thick H II region. We obtain a = —0.59 for 
4C50.55, in agreement with the value given by Mantovani 
et al. (1982) obtained from integrated flux densities com- 
piled from the literature between 178 MHz to 10.7 GHz. 
Spectral indices generally show a relatively flat spectrum 
for the SNR and within the errors of the measurements 
we do not find any variations in the spectral index across 
the remnant, except for the central region (s4), where we 
observe a steeper spectrum. However, in the central re- 
gion the surface brightness is very low, giving rise to high 
uncertainty. 

The boxes el, e2, e3 ,e4 are positioned to sample the 
spectral index of the prominences (see Fig. 3 and Fig. 
5). The general trend is that spectral indices towards the 
prominences are similar to those seen in the main body of 
the SNR's shell, albeit with larger errors arising from the 
weakness of the prominence emission, except for the re- 
gion e2 which indicates a thermal spectrum. We interpret 
the spectral index measurements as supporting evidence 
that three of the prominences indeed are non-thermal and 
related to the SNR, as suggested by the overall SNR mor- 
phology. The region e2, however, seems to be an unrelated 
source of thermal origin, and indeed this prominence is 
unpolarized (see Mantovani et al. (1991) and Section 2.5 
below) . For further evidence of the association of the re- 
maining three prominences to the SNR we must appeal to 
the polarization measurements discussed below. 

The lower-right panel of Figure 4 shows a TT-plot us- 
ing the region indicated by dotted lines in Figure 3. From 
this region we derive a global spectral index of the SNR 
as a ~ —0.43. This calculated spectrum is comparable to 
the value a = —0.42 given by Mantovani et al. (1991), but 
still flatter than the typical value of a = —0.5 for an SNR 
expanding adiabatically. This indicates that the SNR is 
most likely already in the radiative phase. 

2.5. Polarization characteristics of the SNR 

Noise corrected polarized intensity (PI) maps are calcu- 
lated as PI = ^[/2 + Q2 _ (i_2 cr)2, where cr = 30 mK Tg 
is the rms noise in each map. Polarization angle maps have 
been derived, as '0a = 5 arctan ^, from the four bands of 
the DRAO telescopes, each of width 7.5 MHz and centered 
at 1406.65 MHz (band A), 1414.15 MHz (B), 1426.65 MHz 
(C) and 1434.15 MHz (D). These maps have been used to 
calculate the rotation measure across the remnant. In Fig- 



ure 6 we show the polarized intensity and rotation measure 
maps of CTB104A. We observe polarization from almost 
all parts of the remnant. A distinct difference between the 
total intensity and the polarized intensity images is imme- 
diately evident: there are more small-scale structures in 
the polarized intensity image than in the total intensity 
image. This is the result of high rotation measure differ- 
ences on small scales within the remnant, which we assign 
to turbulence and/or localized shock compression. 

The northern part of the remnant is particularly aff'ected 
by depolarization; a distinct region of low polarized inten- 
sity can be seen in this area. Below this area, we ob- 
serve patchy polarization structures. In the polarization 
image it is also possible to see filamentary regions where 
PI is zero. Rather than signalling the absence of polariza- 
tion, these structures are most likely due to beam depolar- 
ization, arising from a fluctuating magneto-ionic medium 
within the beam. 

The percent polarization varies from 5 to 20% over the 
remnant. The extension at the western edge of the SNR 
is weakly polarized. There is polarization towards all 
prominences, except the one in the east. Given the mea- 
sured spectral indices of the prominences (consistent with 
a non-thermal spectrum and also similar to the indices in 
the SNR shell) along with their observed polarization, we 
conclude that with the exception of the eastern one, the 
prominences are associated with CTB104A. 

Towards CTB104A there is smooth detectable total in- 
tensity emission, and yet the detected polarization is domi- 
nated by small-scale structures. Bandwidth depolarization 
can not be the cause, because the four individual bands 
allow us to calculate RM values as high as ~6000 rad m"^. 
The spatial filtering of the smooth structures beyond the 
limited scale of the interferometer affects both the total 
intensity and the polarization channels in the same way. 
However, polarized intensity is additionally affected by the 
foreground Faraday rotation on various scales, even if the 
intrinsic polarization is smooth. 

Clegg et al. (1992) observed RMs from 54 extragalactic 
sources and reported an anomalous RM region, > 1000 
rad mr"^, towards {(., h) ~ (92°, 0°). The angular extent of 
this region is coarsely constrained by point source probes 
to be between 1° and 8°, indicating enhanced RMs behind 
the Cygnus region. By looking at the RM differences of 
the double-lobe sources they conclude that RM differences 
are dominated by sinall-scale structures. The prevalence 
of small-scale variability is obvious in Figure 6, but the 
RM structure also has a systematic component — a grad- 
ual variation of the RMs in the direction of south-east to 
north-west, which is very likely related to the structure 
seen by Clegg et al. (1992). The direction of maximal gra- 
dient in RM lies between 30 ± 10 degrees north of west, 
although this is not easily seen in a gray scale image. 

In order to exhibit this gradient more effectively we 
project the RM distribution on to an axis 30° north of 
west and show the resulting variation of RM along this 
axis in Figure 7. The RM measurements show a system- 
atic increase from negative to positive with respect to the 
position along this axis. We identify this 30° axis as the 
magnetic field axis in the remnant. The observed orien- 
tation of the magnetic field axis supports the idea that 
the shell configuration in CTB104A arises from expansion 
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in a uniform magnetic field. Indeed, this axis, within the 
precision to which it can be ascertained, is the only one 
that can validate this model for CTB104A. The RM gra- 
dient indicates that the magnetic field is tangential across 
the remnant, which is an accepted field configuration for 
relatively old SNRs. 

This RM variation can not be attributed to a fore- 
ground Faraday screen, since the required mean electron 
density would be ~0.2 cm~^ (for an assumed magnetic 
field strength of 1 //G over a pathlength of 1.5 kpc) which 
is extremely high by typical Galactic standards and should 
be observable in emission. Further, a local structure capa- 
ble of producing the RM variation would have to have the 
exact shape and location as GTB104A, which is unlikely. 
This magnetic structure is thus intrinsic to CTB104A or 
to its vicinity. 

3. NEUTRAL MATERIAL NEAR CTB104A 

To search for signatures of interaction of CTB104A 
with the surrounding neutral material we now examine 
the CGPS H I data. The final H I mosaic consists of 9 
CGPS fields, which we have convolved to 2' resolution. 
We have integrated the maps at 2 km s^^ intervals and 
have removed the large scale background emission from 
each channel using the BGF algorithm by Sofue & Reich 
(1979). These maps are plotted in Figures 8 and 9. 

Starting from negative velocities we see emission sur- 
rounding the remnant, but the most important concen- 
tration of H I is around —2 to —10 km s~^, forming a 
D-shaped structure surrounding the SNR. This structure 
slowly fades out with increasing negative velocities. The 
spatial coincidence of the H I with the continuum emission 
is more prominent between —6 and —8 km s~^. 

Between 10 and 8 km s~^ there is an irregular patch of 
H I emission centered on the SNR, which we tentatively 
interpret as the cap of the H I bubble moving away from 
us. The approaching cap is not visible, probably due to 
confusion with material from the local spiral arm. 

As a modifier to the overall shell morphology of the neu- 
tral material surrounding CTB104A, there are clear breaks 
in the shell that correspond well to the positions of two of 
the prominences. Particularly striking is a large break oc- 
curring near the location of the northwest prominence at 
LSR velocities between ~ —14 and —20 km s~^. There 
is a plume of neutral material (at £ ^ 92?6, b ^ G?6) 
which appears to extend beyond, and along a similar axis 
to, the prominence seen at 1420 MHz. In this same ve- 
locity range, the neutral material is evacuated from, and 
appears to encompass, the region into which the south- 
ern prominence extends. Towards the remaining, smaller, 
northeastern prominence associated with the SNR we see 
no clear signs of interaction with the surrounding neutral 
material. However, at the location of this prominence, the 
upper radio shell is broken in two; this suggests that there 
is a similar escape phenomenon occurring here, but that 
the complexity of the H I emission on small scales pre- 
cludes direct identification of interaction of the SNR with 
its surroundings. We thus interpret the prominences as 
hot gas from the interior of the SNR escaping into less 
dense regions, as suggested by Landecker et al. (1985). 

We identify the surrounding H I structure at ~ — 6 
km s~^ as defining the rest velocity (relative to the LSR) 



of CTB104A, which places the remnant at a distance of 
~ 1.5 ± 0.2 kpc, assuming a standard flat Galactic rota- 
tion curve {Rq = 8.5 kpc, Yq = 220 km s~^). At this 
distance the linear size of the remnant is ~35 pc. Note 
that, at lower negative velocities the structure surround- 
ing the northern part of the SNR is more prominent while 
the southern part is best visible at higher negative veloci- 
ties. This indicates that the northern part of the bubble is 
moving away from us while the southern part is approach- 
ing us. Additionally, the radial velocity of the cap at '^lO 
km s~^ indicates an expansion velocity of about 16 km 
for the bubble. 

We also examined our CO-data for a possible correla- 
tion with the SNR. But we could not find any structure 
which might be related. Further observations of the sur- 
roundings and the interior of the SNR are necessary, since 
our measurements only covered the shell of the SNR and 
the eastern H II region. 

4. THE DETECTED H II REGION G94.48-0.3 

Towards £ = 94? 48 and 5 = 0?3 one of the prominences 
is visible in both 408 and 1420 MHz images. In contrast 
to the other three structures, we observe no polarization 
towards this extension of the remnant. Due to the low in- 
tensity values and flat spectral index of the SNR, it is not 
possible to determine whether this prominence is thermal 
or non-thermal. However, examination of infrared data 
shows that this structure is seen as a prominent half-shell 
in 60 /im IRAS emission (see Figure 10). This shell is 
also identified in our CO line data. We have integrated 
the CO emission over the velocity interval —55 km 
to —35 km s~^, over which the CO morphology is simi- 
lar to that of the shell seen at 60 /im (see Fig. 10). The 
coincidence of molecular emission, infrared emission and 
radio emission with a broadly thermal spectrum strongly 
suggests that the shell is an H II region. Moreover, the 
observed LSR velocity of the CO emission indicates that 
the shell has a distance of ~4 kpc. Therefore this previ- 
ously unknown H II region is located behind the supernova 
remnant, which we expected since the polarization images 
reveal no additional depolarization of the supernova rem- 
nant in the direction of the H II region. 

5. DISCUSSION 

CTB104A is striking due to its shape, but more impor- 
tantly due to its distinct rotation measure gradient. Such 
a gradient of RM, with a clear symmetry on both sides of 
the SNR, signifies a tangential magnetic field distribution 
in the remnant. The shock of the explosion compresses 
and deforms the magnetic field and thus the remnant car- 
ries the signature of the local magnetic field. In order to 
interpret this RM gradient in terms of the expansion of the 
SNR in a uniform magnetic field, we need to have informa- 
tion about the large scale magnetic field structure towards 
the SNR. The current method to determine the direction 
of the magnetic field in the Galaxy is the use of pulsar 
rotation measures (see Han et al. 1999), because they are 
relatively easy to measure and in most cases distances of 
the pulsars are known. A plot showing the magnetic field 
direction with respect to the spiral arms of the Galaxy is 
given in Figure 11. According to this figure, the direction 
of the ambient magnetic field is pointing away from us and 



CTB104A: Magnetic Field & Environment 



5 



directed towards lower longitudes. 

Assuming CTB104A is expanding in this uniform large- 
scale magnetic field configuration, it should form a bubble 
by stretching the magnetic field lines as outlined in Fig- 
ure 12. Such a magnetic bubble would result in a decreas- 
ing RM with respect to Galactic longitude (left panel in 
Figure 12). However the observed RM gradient decreases 
in the opposite direction, namely, the RM is increasing 
with longitude. This implies an ambient magnetic field 
direction opposite to the overall Galactic magnetic field. 
Thus the orientation of the magnetic field in this region is 
not parallel to the galactic plane, nor in agreement with 
the "large-scale" Galactic field, but is consistent with the 
RM anomaly observed by Clegg et al. (1992). The over- 
all size of this region is not well constrained by the patchy 
point-source coverage, but it is at least as big as CTB104A 
and may be as large as ~ 8 degrees; such a large scale 
feature must be taken into account when modeling the 
electron distribution of the Galaxy. 

We can make use of the observed variation of the RM 
to estimate the electron density and the magnetic field 
strength in the shell of the SNR. Based on the flat spec- 
trum, we can assume that the SNR has already entered 
the radiative phase. Woltjer (1972) defined the beginning 
of the radiative phase as the time where half of the explo- 
sion energy is lost by radiation. The radius i?rad of the 
SNR at this time depends on its explosion energy Eq and 
the ambient density no before explosion: 

i?rad = 23.1 E^J^' n-'"-'. (1) 
Assuming a typical explosion energy of 10^^ erg, we get a 
lower limit for the ambient density no — 2.0. The com- 
pression of the material 8 in the shell gives an electron 
spectrum, for Fermi acceleration, with 7 defined by 

.^2-^. ,2, 

The spectral index a = —0.43 found in Section 2.4, then 
gives us a (;onipression ratio of 5 = 4.8, which is slightly 
higher than the expected compression of 4.0 for an adi- 
abatically expanding SNR. This results in a lower limit 
to the electron density of nc = 9.6 cm~^, assuming the 
material in the shell is fully ionized. 

The compression of (5 = 4.8 and the linear radius of 
the SNR R = 17.5 pc give the width of the compression 
zone as 1.3 pc. This corresponds to a shell thickness of 
only ^ 3', although the total intensity map gives the im- 
pression of a thicker shell especially towards the north. It 
is not surprising that the thin shell is not noticeable in 
the radio map, because of projection effects. It should be 
remembered that the remnant has not a perfect circular 
shape. It is in an advanced stage of evolution and the am- 
bient medium structure is quite inhomogeneous. Thus the 
expanding shell of the SNR should be locally very thin, 
even though the radio image does not reveal it. 

The rotation measure RM is given by; 

RM = 0.81 ne B|| L (3) 
where is the component of the magnetic field parallel 
to the line of sight and L is the path length through the 
compression zone along the line of sight. For a spheri- 
cal geometry, L can be expressed in terms of the impact 
parameter b as 

L = ^/R^ - 62 - - 62 (4) 



where r is the inner radius of the compressed region. The 
tangential magnetic field configuration then leads to 



So 



(i? + r)/2' 



(5) 



where Bq is the constant magnetic field within the com- 
pression zone. Inserting Eqns 4 and 5 into 3 gives 

RMib) = 1.62 ne Bo , ( \J R^ - \? - y/r"^ - b'^ ) . 

(R + r) \ J 

(6) 

The rotation measure gradient at the center of the SNR 
with respect to the impact parameter, 6, would then be 



dRM{0) 
db 



1.62 ne Bo 



l-r/R 
1 + r/R' 



(7) 



The observed RM distribution, as displayed in Figure 7, 
indicates that the gradient over the inner part of the rem- 
nant is almost constant. We therefore can extract the gra- 
dient given in Eqn 7 by fitting a straight line to the RM 
distribution. The resultant fit gives a slope of 1.30 ± 0.06, 
also plotted in Figure 7. Thus using the compression factor 
S = 4.8, corresponding to r/R = 0.93, and the lower limit 
for the electron density = 9.6 cm~"^ obtained above, 
we calculate the constant magnetic field in the compressed 
region as ~ 2.3 /iG. This would imply an ambient mag- 
netic field strength of about 0.5 jiG. 

6. SUMMARY 

We have presented new, high resolution H I and con- 
tinuum images of CTB104A at 408 and 1420 MHz, aug- 
mented by CO spectral line maps and HIRES IRAS data at 
aOjim.. Despite the superficial complexity of CTB104A, a 
simple physical interpretation of the dominant processes at 
work in the SNR emerges from the CGPS data. We iden- 
tify the orientation of the local magnetic field, the physical 
origin of the shell emission and the overall morphology of 
the surrounding neutral shell that encloses the SNR. The 
remnant as a whole is indeed complex, however, as re- 
vealed by the signatures of localized interaction between 
the remnant and the surrounding neutral medium. 

Polarization and rotation measure structures are dom- 
inated by small scale variations, induced by turbulence 
within the remnant. However, there is a distinct, or- 
dered, rotation measure gradient from negative values in 
the southeast to positive ones in the northwest. The di- 
rection of the rotation measure gradient defines the axis of 
the magnetic field around the remnant, lying 30° north of 
west. This magnetic field is local to the remnant and gen- 
erally counter to the global magnetic field of the Galaxy 
in this area. This magnetic field structure lies inside a 
large RM anomaly inferred from rotation measure data by 
Clegg et al. (1992). The magnetic field strength within 
CTB104A is Bo ~ 2.3 /xG. 

Non-thermal radio continuum emission arises from two 
shell edges at the expected locations of compression for 
this magnetic field configuration. A broken shell of neutral 
material surrounds the SNR. One of the continuum shell 
edges, to the northeast, is broken in two. At the location 
of the break, the smaller of CTB104A's three prominences 
extends from the shell. We find no obvious disruption, or 
break, in the surrounding neutral shell at the location of 
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this smaller prominence, probably due both to the smaller 
size of the prominence and the overall complexity of the 
surrounding H I emission. Towards the two larger promi- 
nences, however, the neutral shell is clearly broken, having 
been presumably less dense at these points prior to disrup- 
tion. The morphology of the neutral medium, as traced 
by H I, is consistent with it being pushed by the escaping 
ionized material in these regions. 

In summary, the three prominences around the remnant 
reveal that CTB104A is probably the best example of the 
leakage of blasts of ionized material into a low density envi- 
ronment. One extension to the east is identified as a back- 
ground H II region associated with a massive molecular 
shell. CTB104A is a relatively old supernova remnant, as 
deduced from its global spectral index, expanding isother- 



mally within a uniform magnetic field and inside an asso- 
ciated H I shell. The distance 1.5 kpc inferred from the 
radial velocity of the H I, Vlsr — —6 km s~^, translates 
to a mean linear diameter of ~ 35 pc. 
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CTB104A: Magnetic Field & Environment 



Fk;. 1. — 1420 MHz contour plot (top) and gray scale image of CTB104A with a resolution of 65" x 49" (EW x NS). Contour levels start 
from 4.8 K Tb and run in steps of 0.5 K Tb. The diffuse structure with a bright core in the lower left corner is the H II region S124. 



Fig. 2. — 408 MHz contour plot (top) and gray scale image of CTB104A with a resolution of 3'.7 x 2.'8 (EW x NS). Contour levels start 
from 45 K Tb and run in steps of 10 K Tb. The diffuse structure with a bright core in the lower left corner is the H II region S124 



Fig. 3. — Regions used for the TT-plot analysis are displayed on a gray scale image of the region at 1420 MHz on a grid of 1 arc minute. 
The circles show the regions around the point sources which were excluded. The box with dotted lines shows the region used for an overall 
spectral index calculation. The dashed circle,centered at £ = 93? 7 ad 6 = — 0?3, shows the approximate circular part of the remnant and has 
angular dimensions of 1?2 X 1?2. The arrows indicate the positions of the prominences extending away from the remnants boundary 



Fig. 4. — Results of the TT-plot analysis for the selected regions, numbered from 1 to 6. The spectrum relevant for S 124 and 4C50.55 are 
given in the lower panel. The lower-right panel shows the spectra of the SNR obtained from the region plotted with dotted lines in Figure 3 



Fig. 5. — Results of the TT-plot analysis towards the four extensions of the SNR 



Fig. 6. — Polarized intensity image (top) of CTB104A at 1420 MHz with overlaid total intensity contours. Contours start at 7.0 K Tb and 
run in steps of 0.5 K Tb. Lower panel shows the rotation measure image of the same region. White contour is at —100 rad m~^ and black 
contour is at 100 rad m~^ 



Fig. 7. — Rotation measure distribution as a function of the angular distance from the center of CTB104A, (£,b)= 93?7, — 0?3, projected 
on to the magnetic field axis, which maJies an angle of 30° with respect to the Galactic plane. The straight line is the fit to the RM gradient 
as explained in the text 



Fig. 8. — Neutral hydrogen data at 2' resolution integrated over 2 km s ^ intervals, between 10 and —6 km s ^. The overlaid contours are 
from the 1420 MHz continuum image 



Fig. 9. — Same as Figure 8 but for the velocity interval —8 to —24 km s 



Fig. 10. — CO emission over the velocity interval between —55 and —35 km s ^ (left) and infrared 60 /im (right) images towards CTB104A. 
White contours show the total intensity emission at 1420 MHz 



Fig. 11. — Direction of the magnetic field with respect to the spiral arms in the Galaxy, obtained from the pulsar rotation measure data, 
adopted from Han et al. (1999). The approximate position of CTB104A is marked by a cross 



Fig. 12. — Plot of the model for an SNR expanding in a regular magnetic field. The magnetic field direction and its parallel component 
given by the pulsar rotation measure data arc shown on the left. Also shown are the stretch of the magnetic field lines due to the SNR 
expansion and the relative RM strength and direction with respect to the given magnetic field configuration. The sketch on the right shows 
the same configuration but the direction of the magnetic field is reversed 
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